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I INTRODUCTION 


The  purpose  of  the  research  contract  entitled  "Measurement  of  Sea 
Scatter  and  Buoy  Tracks  at  Long  Ranges  by  High-Resolution  OTH-B  Radar" 
is  to  determine  the  feasibility  of  using  a skywave  HF  radar  to  remotely 
measure  the  ocean  wave  spectrum,  the  oceanic  wind  field,  and  the  ocean 
currents  in  the  North  Pacific  Ocean.  The  work  was  performed  over  a two- 
year  period  extending  from  1 December  1973  through  31  December  1975. 

Our  work  involved  the  cooperation,  support,  and  coordination  of  the 
scientists  at  Scripps  Institution  of  Oceanography  (SIO);  Stanford 
University  ( SU ) ; and  Texas  A & M University  (TAMU). 

This  report  consists  of  a brief  outline  of  our  work  (Sections  II 
through  IV),  and  copies  of  papers  submitted  or  accepted  for  publication 
(Appendices  A,  B,  and  C). 
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II  MEASUREMENT  OF  SEA- SCATTER  AND  BUOY  TRACKS 
BY  HIGH -RESOLUTION  OTH-B  RADAR 


SRI's  participation  during  the  first  year  of  the  NORPAX  effort  was 
coordinated  with  and  in  si^pport  of  the  more  basic  oceanographic  research 
by  the  scientists  at  SIO,  SU,  and  TAMU.  This  effort  was  divided  into 
two  tasks.  The  first  task  was  the  recording  and  spectral  analysis  of 
sea  backscatter  data  over  predesignated  areas  in  the  North  Pacific  Ocean. 
Doppler  spectra  were  given  to  the  scientists  at  SIO  for  analysis  of  wind 
velocity.  The  second  task  was  the  recording  and  analysis  of  backscatter 
range-azimuth  tracking  data  to  determine  the  position  of  a network  of 
drifting  buoys,  each  instrumented  with  an  HF  repeater. 

A.  Tracking  Drifters 

The  tracking  of  drifting  buoys  was  begun  during  the  NORPAX  Pole 
Experiment  and  was  successfully  completed  with  an  rms  position  error 
of  20  km.  Six  drifting  buoys,  three  drogued  at  30  m with  8.5-m-diameter 
parachutes  and  three  undrogued,  were  launched  January  29,  1974.  Two  of 
the  drifters  (one  drogued  and  one  undrogued)  malfunctioned  immediately ; 
two  of  the  drifters  (both  undrogued)  failed  early  in  the  experiment  or 
rapidly  drifted  out  of  the  radar  coverage;  one  drifter  (drogued)  was 
tracked  for  one  month,  after  which  it  failed;  and  one  drifter  (drogued) 
was  tracked  for  four  months.  The  results  of  this  effort  are  found  in 
a report  by  Marshall  and  Barnum.1*  The  tracking  Information  was  given 


References  are  listed  at  the  end  of  the  main  text.  References  for 
the  appendices  are  listed  at  the  end  of  each  appendix. 
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to  the  scientists  at  SIO  and  TAMU  for  interpretation.  The  results  and 
interpretation  of  the  long-term  trajectories  are  presented  by  Kirwan 
and  McNally. 

B.  Sea  Backscatter  Data 

Recording  and  processing  of  the  backscatter  measurements  were  also 
begun  during  the  NORPAX  Pole  Experiment.  The  results  of  this  effort  are 
also  found  in  Marshall  and  Barmim/  Both  wide-area  scans  in  the  vicinity 
K.V.  Flip  and  longer  dwells  over  specific  areas  of  the  North  Pacific 
were  made.  The  wide-area  scans  were  processed  for  12.8  s of  coherent 
integration  time  for  21  range  lines  (1.5  km  per  range  line),  and  the 
longer  dwells  were  processed  for  25.6  and  51.2  s of  coherent  integration 
time  over  8 range  lines.  The  longer  dwells  were  taken  each  day  between 
2J00  and  2400  GMT  when  the  HF  repeater  aboa'd  R.V.  Flip  was  turned  on. 

The  longe r-dwe 1 1 sea  backscatter  data  taken  at  this  time  were  not  analyzed 
because  of  multipath.  Better  spectra  could  have  been  obtained  by 
operating  at  a different  time  of  day.  For  the  wide-area  scans  taken 
during  the  Pole  Experiment  and  during  a March  1974  storm  passage  west  of 

O 

150  W,  the  sea  backscatter  Doppler  spectra  were  obtained  for  analysis  and 
correlation  with  the  observed  wind  speed  and  direction  obtained  from 
R.V.  Flip  and  other  ship  reports.  These  data  were  given  to  the  scientists 
at  SIO  for  analysis  and  comparison  with  the  ground  truth.  The  results  of 
the  correlations  between  wind  speed  and  direction  have  been  presented  by 
Stewart  and  Barnum. 
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Ill  REMOTE  SENSING  OF  SEA  SURFACE  MOVEMENT 
AT  LONG  RANGES  BY  IONOSPHERICALLY  PROPAGATED 
BACKSCATTER  WITH  HIGH  RESOLUTION 

The  second  year  in  the  NORPAX  effort  was  a coordinated  undertaking 
between  the  scientists  at  SIO  and  SU.  The  purpose  of  this  research  was 
to  record,  spectral-analyze,  and  correlate  the  backscatter  measurements 
with  observed  values  of  surface  winds  and  surface  waves  obtained  from 
ships  of  opportunity.  As  part  of  this  effort,  basic  questions  about 
the  optimal  spectral  processing  necessary  to  interpret  the  sea  back- 
scatter data  were  investigated. 

The  results  are  presented  in  three  papers  (Refs.  4-6).  Wind  speed 
estimates  were  improved  by  normalizing  the  data  by  the  inverse  of  the 
square  root  of  the  frequency.  Accuracies  of  2.5  m/s  were  obtained. 

The  results  of  the  wind-speed  correlation,  the  applicability  of  the 
technique,  and  a discussion  of  the  potential  errors  are  discussed  by 
Maresca  and  BarnunT  in  a paper  accepted  for  publication  in  1977,  further 
extending  the  discussion  in  Ref.  3. 

The  ability  to  map  the  entire  wind  direction  field  is  discussed  by 
Barnum,  Maresca,  and  SerebrenyL  in  a paper  accepted  for  publication  in 
1977.  This  paper  is  included  in  our  final  report  as  it  is  directly 
pertinent  to  NORPAX  work  even  though  the  sea  backscatter  data  were  taken 
prior  to  our  contract.  The  high-resolution  measurements  in  the  vicinity 
and  across  the  atmospheric  front  are  especially  noteworthy. 

Ionospher ica 1 ly  clean  Doppler  spectra  were  obtained  for  coherent 
integration  times  of  102  s.  Johnstone's7  second-order  radio-wave 
scattering  model  for  bistatic  radars  was  used  in  a trial-and-error 
calculation  to  determine  the  ocean  wave  spectrum  and  is  discussed  by 
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Maresca,  Johnstone,  and  Bamum  in  a paper  to  be  submitted  to  the  Journal 
of  Geophysical  Research.  Ocean  wave  spectra  were  not  available  for 
comparison,  but  hourly  observation  of  the  wave  height  by  two  Navy  ships 
in  the  area  agreed  within  107.  of  the  significant  wave  height  calculated 
from  the  radar-derived  wave  spectra. 
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IV  PRINCIPAL  RESULTS 


High-frequency  skywave  radar  radio  wave  measurements  of  wave  height 
and  wind  speed  are  consistent  with  theory. 

The  ability  to  scan  large  areas  allows  the  derivation  of  the  entire 
wave  and  wind  field  for  large  storm  passages.  The  high  resolution  of 
WARF  allows  investigation  of  intense  gradients  across  low-pressure  cells 
and  across  atmospheric  fronts. 

Wind  speed  estimates  were  accurate  to  2.5  m/s  and  wind  direction 
estimates  were  accurate  to  16°.  Wave  height  measurements  derived  from 
the  radar-measured  ocean  wave  spectrum  were  accurate  to  within  10%  of 
the  observed  wave  heights.  Skywave  radar  position  fixes  used  to  measure 
ocean  currents  by  tracking  drifting  buoys  were  accurate  to  within  20  km. 

The  ionosphere  will  support  propagation  sufficiently  to  make  the 
same  type  of  measurements  as  line-of-sight  radar  with  approximately  the 
same  accuracies.  The  important  difference  is  the  occasional  degradation 
of  skywave  data  by  unstable  propagation. 

For  a limited  data  set,  it  was  found  thac  the  incoherent  averaging 
of  A to  8 Doppler  spectra  obtained  at  the  same  range  line  or  the  incoherent 
averaging  of  l to  2 Doppler  spectra  obtained  over  an  8-range-line  average 
will  produce  a stable  mean;  the  ratio  of  the  Bragg  lines  was  used  as  a 
parameter.  Increasing  the  coherent  integration  time  did  not  significantly 
affect  the  mean  but  did  allow  for  greater  resolution  in  the  Doppler 
spectra . 
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Appendix  A 


HIGH -RESOLUTION  MAPPING  OF  OCEANIC  WIND  FIELDS 
WITH  SKYWAVE  RADAR 


ABSTRACT 

The  direction  of  the  mean  surface  wind  field  in  the  North 
Pacific  Ocean  was  mapped  on  25  and  26  September  1973  over  an 
area  of  3 X 10^  (km)2  by  OTH-B  HF  radar.  A spatial  resolution 
of  60  km  in  range  and  15  km  in  cross  range  were  used  at  points 
spaced  by  150  km  in  range  and  80  km  in  cross  range.  Wind  di- 
rections were  inferred  from  the  upwind/downwind  first-order  Bragg 
ratio  and  the  measure  of  the  maximum  ratio  occurring  for  radial 
winds  at  points  near  each  observation.  Detailed  measurements 
of  the  wind  speed  and  direction  were  used  to  define  the  position 
of  an  atmospheric  cold  front.  Range  resolutions  of  3 km  are 
achieved  with  the  WARF  radar  and  are  essential  to  the  analysis 
of  atmospheric  fronts  and  the  gradient  of  the  wind  stress  across 
intense  cyclones. 


I INTRODUCTION 

Measurements  of  ocean  surface  winds  are  important  to  the  under- 
standing of  air/sea  interaction.  Along  with  desert  and  polar  areas, 
the  oceans  comprise  vast  data-sparse  regions.  The  violent  and  de- 
structive tropical  cyclones  are  spawned  and  grow  to  great  dimensions 
and  intensities  in  such  areas.  On  these  occasions,  ocean  surface 
measurements  are  esp  dally  desirable. 

Inferences  of  ocean  surface  waves  and  winds  have  been  achieved 
remotely  with  the  use  of  high-frequency  (HF)  radar  [l]-[8].  Radar 
echoes  from  the  sea  at  decameter  wavelengths  are  caused  by  Bragg 
resonance  of  the  radio  waves  with  individual  ocean  waves  and  groups 
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of  interacting  waves  [3], [6].  First-order  Bragg  resonance  with  ocean 
waves  of  one  half  the  radio  wavelength  produce  the  strongest  scatter. 

With  skywave  [2]  (ionospheric)  propagation  at  6 to  30  MHz  those  waves 
have  4.0-to-1.8-s  periods,  respectively,  and  are  tightly  coupled  to  the 
local  winds  especially  for  frequencies  above  10  MHz.  Higher-order  scatter 
is  produced  by  ocean  waves  of  all  lengths  greater  than  half  the  radio 
wavelength,  and  the  most  significant  returns  are  found  at  Doppler  fre- 
quencies near  the  first-order  Bragg  return  [5]. 

Long  and  Trizna  [1]  were  first  to  publish  radar  maps  of  ocean  wind 
direction  obtained  at  long  ranges  with  skywave  propagation  [2],  Barrick 
e_t  £l . [3]  inferred  significant  wave  height  from  surface  wave  radav  echoes 
at  close  ranges.  Ahearn  e£  al^.  [4]  mapped  ocean  wind  directions  with  the 
method  in  [1],  and  also  attempted  to  infer  ocean  wind  speeds  from  the 
ocean  scatter  magnitude  near  zero  Doppler.  Tyler  £t  aK  [5]  inferred 
the  directional  spectrum  of  sea  waves  near  Wake  Island  with  a bistatic 
surface-wave  radar.  Johnstone  [6]  published  a detailed  treatment  of  the 
first-  and  second-order  theory  of  HF  scatter  from  the  sea,  and  found 
agreement  in  the  measure  of  wind  speed  with  surface-wave  sea-scatter 
recordings  near  Wake  Island  during  the  same  experiment  [5]. 

A new  method  for  inferring  wind  fields  at  long  ranges  by  skywave 
radar  was  presented  by  Stewart  and  Barnum  17].  Wind  direction  was 
inferred  from  the  first-order  Bragg  scatter  ratio  from  approaching  and 
receding  waves,  using  a new  model  for  the  ocean  wave  directional  spectrum 
that  is  a function  both  of  wind  speed  and  radio  frequency.  The  wind 
speed  was  inferred  from  the  Doppler  spectrum  width  centered  on  the 
strongest  Bragg  return  at  a level  10  dB  lower  than  the  Bragg  peak. 

Maresca  and  Barnum  [8]  have  improved  this  measure  of  wind  speed  by 
accounting  for  the  effect  of  radar  frequency  and  have  discussed  the 
theoretical  limitations.  Reliable  measures  of  ocean  surface  winds 
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accompanied  the  data  discussed  in  [7]  and  [8],  from  which  we  now  deduce 
accuracies  of  ±16°  in  wind  direction  and  ±2.4  m/s  in  wind  speed  within 
one  standard  deviation  relative  to  the  measures  i_n  situ. 

Skywave  radar  coverage  far  exceeds  that  of  HF  surface  wave.  HF 
radar  analysis  of  large  weather  disturbances  such  as  storm  systems  would 
be  possible  only  with  the  larger  coverage.  By  comparison,  the  use  of 
surface-wave  radar  has  contributed  most  extensively  to  our  understanding 
of  the  complete  Doppler  frequency  spectrum  of  sea  scatter  [3), [5], [6]. 
Since  the  lower  atmosphere  has  negligible  effect  on  the  coherence  of 
the  propagating  wave,  and  only  a single  mode  (the  surface  wave)  is 
present,  very  long  coherent  data  integrations  are  thereby  permitted 
(e.g.,  tens  of  minutes).  Hence,  the  Doppler  resolution  can  be  extremely 
high.  Such  integration  followed  by  several  incoherent  spectrum  averages 
yields  a smooth  estimate  of  the  mean  spectrum  amplitudes  of  the  scatter. 
Essentially,  two  potential  problems  exist  with  the  use  of  ionospheric 
propagation  for  the  purpose  of  sea  scatter  analysis:  (1)  a coherence 

time  limit  caused  presumably  by  the  motion  of  fine-scale  irregularities; 
and  (2)  the  simultaneous  reception  of  more  than  two  propagating  modes 
(multipath)  from  different  areas  of  the  sea.  The  coherence  time  allowed 
by  the  ionosphere  will  generally  limit  the  Doppler  resolution  to  between 
0.1  and  0.01  Hz  [8], [9]  at  middle  latitudes,  and  can  be  expected  to  be 
poorer  at  more  northerly  latitudes  [9],  The  highest  ■'herence  is  found 
for  E-layer  reflection,  because  it  affords  the  least  penetration  of  the 
ionosphere.  Unfortunately,  this  low  layer  prohibits  propagation  beyond 
about  2000  km  range  via  a single  hop,  and  it  is  generally  absent  between 
sunset  and  sunrise.  Multipath  propagation  is  evident  in  a large  percen- 
tage of  data  recorded  [8]  and  is  considerably  more  prevalent  at  ranges 
exceeding  3000  km  via  single-hop  F-layer  propagation  (7).  Multipath 
will  always  obscure  some  portion  of  the  sea  scatter  Doppler  spectrum, 
when  viewed  as  a whole,  but  will  not  necessarily  limit  inference  of 
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wind  direction  or  speed  from  the  scatter  located  near  the  two  first-order 
Bragg  returns  [8]. 

This  paper  illustrates  the  results  of  a two-day  skywave  radar 

scatter  experiment  over  the  Pacific  Ocean.  The  direction  of  the  ocean 

6 2 

wind  field  was  mapped  to  high  resolution  over  3 X 10  km  and  the  wind 
velocity  was  detailed  across  a cold  front.  It  will  be  shown  that 
significant  changes  in  wind  direction  occur  over  space  scales  on  the 
order  of  20  km  or  less,  such  that  high  radar  resolution  was  indispensible . 
Ionospheric  effects  on  these  data  are  discussed. 


II  THEORY 

The  method  for  determination  of  wind  speed  from  the  second-order- 
scatter  contribution  to  the  Doppler  spectra  has  been  discussed  elsewhere 
[7], [8]. 

The  determination  of  wind  direction  for  the  September  1973  experiment 
is  analogous  to  that  of  Long  and  Trizna  [1],  To  infer  the  wind  direction 
we  assumed  that  the  2-to-3-second  waves  responded  within  a short  time 
(approximately  one  hour)  to  the  local  winds,  and  that  the  measure  really 
represents  the  mean  wind  over  the  same  period  of  time. 

The  sea-wave  energy  directional  spectrum,  F(9),  for  (d | £ r/2  was 
assumed  to  be  somewhere  between  semi-isotropic  and  cosine  squared: 

2 

F(9)  = A + B cos  9 (1) 

where  A and  B are  arbitrary  constants,  and  9 is  the  angle  between  the 
wind  and  wave  directions.  It  was  also  found  empirically  [1]  that 


F (t 


- 0)  = 


F2(tt/2) 

F(9) 


(2) 
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We  define  the  ratio  of  first-order  Bragg  line  amplitudes  in  the  approach- 
ing and  receding  directions  by  X(8),  where 


X (9 ) 


A F (9  ) 

" F (jt  - 9) 


(3) 


Nearest  each  data  scan  we  measure  the  maximum  X(8)  at  8 * 0 or  the  inverse 
of  the  minimum  X at  9 = ~ , with 
X(~/2)  = 1.  We  then  find  that 


of  the  minimum  X at  9 = ",  with  either  one  called  X , and  we  define 

m 


X (6 ) 


fu  + K cos2  8]  . . . 0 £ 1 9 1 * tt/2 

2 -2 

[1  + K cos  6]  . . .^t/2  £ |e  | £ 7T 


where 


k - A TT 


m 


(4) 


(5) 


Equation  (4)  can  be  inverted  to  express  8 as  a function  of  X(8). 

In  [1]  and  [4]  we  find  that  X was  assumed  to  be  constant  for  all 

m 

data  recorded,  regardless  of  the  operating  frequency  or  wind  speed.  In 

light  of  the  work  by  Tyler  et  al.  [5]  we  cannot  hold  X constant  for  our 

m 

experiment.  By  assuming  constant  X we  fix  the  sea  wave  directional 

m 

spectrum  for  all  time.  Tyler  et^  al_.  concluded  that  F(9)  is  a function 
of  both  the  ratio  frequency  and  wind  speed,  and  use  a different  func- 
tional representation  for  F(6),  which  is  a modification  of  the  form 
suggested  by  Longuet-Higgins , Cartwright,  and  Smith  [10].  When  wc 

constrain  the  F(9)  in  [5]  to  match  observed  values  of  X and  X(n/2)  * 1, 

m 

we  calculate  a maximum  of  only  11°  difference  between  9 predicted,  using 

F(6),  by  both  methods.  We  have  used  the  measured  X nearest  to  each 

m 

recorded  spectra,  thereby  updating  Eq.  (4)  with  Eq.  (5)  for  all  data 

reported  herein.  X varied  between  14  and  24  dB,  with  a mean  of  21  dB. 

m 

We  believe  that  any  effects  on  F(8)  caused  by  changes  in  radio  frequency 
and  wind  speed  have  thus  been  partially  accounted  for. 


17 


\ 

III  RADAR  CHARACTERISTICS 

The  Wide  Aperture  Research  Facility  (WARF)  was  designed  originally 
for  high-resolution  DOD-related  applications,  and  has  not  been  described 
to  any  great  detail  in  the  open  literature.  The  reader  is  referred  to 
a paper  by  Skolnik  and  Headrick  [11]  for  an  introduction  to  OTH  radar 
practices,  and  to  a paper  by  Crofc  [2]  for  illustrations  of  data  recorded 
via  skywave  propagation.  Some  details  regarding  the  WARF  radar,  and 
examples  of  swept-frequency  ocean  backscatter  amplitude  were  publi~hed 
by  Barnum  [12].  A more  recent  summary  was  published  in  report  form  by 
Marshall  and  Barnum  [13].  The  WARF  coverage  is  shown  in  Fig.  1. 

Briefly,  WARF  is  a bistatic  system  with  185  km  separation  between 
transmitter  and  receiver  in  central  California.  Four  megawatts  of 
average  effective  radiated  power  (ERP)  are  directed  toward  the  scatter- 
ing area.  This  large-receiving-aperture  (2.55  km)  radar  yields  a radar 
antenna  beamwidth  of  1/3°  at  20  MHz.  A time-delay  resolution  of  20  us 
was  used  to  record  our  data  (where  1 us  is  close  to  150  m of  ground  range 
for  this  skywave  experiment).  Thus,  we  processed  ocean  areas  with  dimen- 
sions of  15  km  (average)  in  cross  range  and  3 km  in  range.  A total  of 
21  range  lines  were  averaged  incoherently  to  smooth  the  scattering 
statistics,  resulting  in  a net  spatial  average  of  15  km  X 60  km. 

Performance  assurance  of  the  WARF  backscatter  system  is  achieved 
through  simultaneous  operation  of  two  additional  systems.  The  first 
system  provides  continuous  surveillance  of  the  HF  operating  frequencies 
available  for  propagation  to  ensure  that  interfering  signals  that  can 
confuse  the  spectral  information  in  the  backscatter  are  not  present  in 
the  swept  bandwidth.  The  second  system  is  a low-power  wide-t requency 
swept-bacL scatter  system  used  to  estimate  the  mode  of  signal  propagation 
to  the  area  of  interest  and  to  specify  usable  operating  frequencies. 
Examples  of  such  data  were  published  by  Croft  [2]  and  Barnum  [12]. 
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IV  SPECTRAL  PROCESSING 


Figure  2 illustrates  the  basic  form  of  data  output  from  the  radar 
processor.  Three-dimensional  displays  of  backscatter  amplitude  versus 
time  delay  (range)  and  Doppler  frequency  are  first  processed  by  means 
of  a double  Fourier  transform  of  the  deramped  SFCW  signals.  The  Doppler 
coverage  extends  from  +2.5  Hz  to  -2.5  Hz  for  an  SFCW  repetition  rate  of 
5 Hz.  The  raw  samples  are  weighted  by  a cosine-squared  window  to  reduce 
frequency  sidelobes.  A coherent  integration  of  12.8  s was  used  exclusive- 
ly with  the  data  discussed  in  this  paper. 


FIGURE  2 EXAMPLE  OF  THE  OCEAN-WAVE  BACKSCATTER 
DOPPLER  SPECTRA.  Bottom  is  spectrum  vs  range 
(time  delay)  at  a single  azimuth.  Top  is  the  average 
of  a'l  range  lines. 
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SPECTRUM  AMPLITUDE 


Incoherent  spectrum  averaging  is  used  to  smooth  the  statistical 
variation  of  sea  scatter  amplitude.  This  variation  was  found  to  resemble 
a Gaussian  process  both  in  time  and  space.  Two  incoherent  time  averages 
and  21  spatial  (range-cell)  averages  were  performed  prior  to  each  spectrum 
output  in  this  experiment  and  should  produce  a very  smooth  estimate  of 
the  representative  scatter  spectrum  1 8 ] . 

Figure  3 illustrates  the  portion  of  the  real-time  processor  output 
used  in  this  experiment.  The/e  data  were  updated  once  per  minute.  Three 
adjacent  antenna  beams,  spaced  by  one-half  degrep,  were  sampled  and 
processed  simultaneously.  A log  of  radar  operating  parameters  was 
recorded  automatically.  The  Bragg  line  ratio  X,  and  spectrum  width  B 
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FIGURE  3 EXAMPLE  OF  THE  REAL-TIME  BACKSCATTER  DISPLAY  USED  FOR 
RECORDING  DATA  FOR  WIND-FIELD  MAPPING.  Each  such  display 
is  updated  once  per  minute  at  a new  ocean  location.  The  center  azimuth 
was  used  to  map  wind  fields. 
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were  also  calculated  automatically  for  each  azimuth,  and  are  labeled 
X and  B,  respectively,  on  Fig.  3. 


V DATA  COLLECTION  AND  INTERPRETATION 

Data  similar  to  Fig.  3 were  recorded  sequentially  by  stepping  every 
2°  In  azimuth  at  a constant  radar  time  delay  and  frequency.  The  frequency 
was  changed  prior  to  each  scan  to  ensure  adequate  propagation  along  each 
scan.  Time  delays  were  stepped  in  1-ms  increments  between  10  and  27  ms 
on  the  first  day  and  between  13  and  25  ms  on  the  second  day  using  F2- 
layer  propagation.  The  longer  ranges  were  sampled  earlier  in  the  dry 
because  multipath  was  least  prevalent  at  those  times. 

Propagation  conditions  were  usually  different  at  tbe  two  ends  of 
an  azimuth  sector,  such  that  more  optimum  time  delays  and/or  frequencies 
could  have  been  chosen  for  different  portions  of  each  scan.  Ionospheric 
multipath  often  appeared  over  portions  of  a scan  at  azimuths  where  the 
skip  distance  was  approached  too  closely.  The  calculation  of  X(6)  for 
wind  direction  was  rarely  affected  by  multipath,  but  not  all  spectra  can 
be  reduced  for  wind  speed  using  B.  The  data  beyond  20  ms  (2000  km)  were 
the  most  often  perturbed. 

Figure  4 illustrates  spectra  from  a portion  of  one  such  azimuth  scan 
on  26  September  1973.  These  data  illustrate  three  results  from  this 
experiment:  the  occasional  observation  of  multipath,  as  shown  at  -6° 
azimuth  (276°  True  from  the  radar);  a double  wind  shift,  noted  by  reversals 
in  the  sign  of  X across  a weather  front;  and  a more  gradual  wind  direction 
change  SW  of  a high  pressure  area.  Multipath  occurs  on  other  spectra  to 
a lesser  degree  (e.g.,  at  -4°)  but  did  not  degrade  the  calculation  of  B. 

All  spectra  from  the  experiment  were  first  reduced  for  the  Bragg 
line  ratio,  X,  and  this  was  plotted  vs.  azimuth,  parametric  in  time 
delay.  Figure  5 is  an  example  of  some  of  these  plots  on  26  September 
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FIGURE  5 


BRAGG-LINE  RATIO,  X(0)  vs  RADAR  AZIMUTH  AT  EIGHT  TIME  DELAYS  ON 
26  SEPTEMBER  1973.  Circulation  to  the  SW  of  a high-pressure  system  was  sampled 
by  the  data  between  13  and  16  ms,  while  a cold  front  was  crossed  by  the  scans 
between  21  and  24  ms. 


1973.  The  plots  on  the  left  of  the  figure  represent  scans  to  the  west 
of  an  anticyclone  and  to  the  east  of  a cold  front.  The  maximum  observed 
values  of  X (in  absolute  dB)  occurred  in  the  trade  wind  region,  since 
the  winds  were  approximately  parallel  to  the  radar  heading.  The  plots 
on  the  right  side  of  the  figure  contain  radar  crossings  of  a weather 
front  that  noticeably  interrupted  the  otherwise  continuous  nature  of 
X versus  azimuth.  The  weather  front  usually  caused  a double  zero  cross- 
ing in  X that  is  seen  at  azimuths  between  -6°  and  2°  in  Fig.  5,  depending 
on  the  time  delay.  These  zero  crossings  essentially  accounted  for  net 
changes  in  wind  direction  exceeding  180°.  Note  that  the  minimum  in  X 
between  these  crossings  advances  in  azimuth  from  -4°  at  21  ms  delay, 
to  between  0°  and  2°  at  24  ms  delay.  This  trend  was  interrupted  by  a 
dip  at  -2°  at  22  ms  with  a lack  of  a zero  crossing.  Since  crossings 
were  also  observed  at  19  ms  and  20  ms  delays,  a reversal  of  wind  direc- 
tion was  therefore  assumed  to  occur  at  -2°  azimuth  at  22  ms. 

The  left/right  ambiguity  in  9 [1],{4]  was  resolved  by  general 
knowledge  of  the  overall  circulation  around  the  high-pressure  cell 
shown  on  the  NWS  map  and  cloud  photograph.  Each  zero  crossing  yields 
a wind  that  is  crosswise  to  the  radar  bearing,  a condition  that  usually 
occurred  across  the  cold  front. 


VI  RADAR  WIND  FIELD  MAPS 

The  radar-derived  wind  field  for  the  first  day,  25  September  1973, 
is  shown  in  Fig.  6,  and  is  compared  with  the  relevant  portion  of  the 
National  Weather  Service  (NWS)  surface  analysis  for  1800  GMT.  As 
indicated  on  the  table  in  Fig.  6,  the  radar  data  were  recorded  between 
the  hours  of  1855  GMT  and  2512  (0112)  GMT.  We  expect  the  data  at  the 
lower  time  delays  (later  in  the  day)  to  disagree  with  the  ship  reports 
taken  some  6 to  7 hours  earlier. 


25 


26 


RADAR  WEATHER  MAP  FOR  25  SEPTEMBER  1973,  COMPARED  TO  THE  NWS  ANALYSIS  OF  1800  GMT 
ON  THE  SAME  DAY.  The  radar  winds  represented  are  for  times  between  1800  and  2400  GMT.  Wind  direction 
is  plotted  by  arrows;  the  feet  of  the  arrows  denote  positions.  Wind  circulation  generally  agrees  well  with  the 
weather  map.  The  cold  front  detected  by  the  radar  does  not  agree,  however. 


The  ESSA-VII  cloud  photograph  taken  near  1800  GMT  is  shown  In 
Fig.  7 for  comparison  to  Fig.  6.  The  cold  front  is  identified  by  a 
well -developed  elongated  bright  band  of  clouds  oriented  NE/SW  extending 
through  40°N,  150°W.  The  cellular-appearing  clouds  between  140°W  and 
130°W  and  near  30°N  are  characteristic  of  clouds  around  the  southwestern 
perimeter  of  an  oceanic  high-pressure  area  ahead  of  a well  developed 
cold  front.  The  center  of  the  high-pressure  area  is  located  somewhere 
in  the  clear  area  in  the  right-hand  edge  of  the  photograph  (Fig.  7). 

The  general  agreement  between  the  radar-derived  wind  fields  and 
the  independent  surface-wind  measurements  on  25  September  was  good. 

The  cold  front  position  shown  in  Fig.  6 was  plotted  by  NWS  approximately 
200  nmi  to  the  northwest  of  the  front  indicated  by  the  radar  data.  The 
radar- inferred  wind  shifts  agree  with  the  wind  shifts  implied  by  the 
satellite  cloud  photograph  displaying  and  locating  the  frontal  cloud 
cover.  The  discontinuous  nature  of  the  wave  field  recorded  by  the  radar 
is  a good  method  for  detailing  the  position  of  the  front  and  the  wind 
field  across  the  front. 

The  radar-derived  wind  field  for  the  second  day,  26  September  1973, 
is  shown  in  Fig.  8,  and  the  corresponding  cloud  photograph  in  Fig.  9. 

The  radar  winds  showed  the  frontal  wind  discontinuity  very  clearly,  as 
well  as  a wind  discontinuity  associated  with  a trough  in  the  easterlies 
west  of  155  °W  and  proximal  to  25 °N  (Fig.  8).  With  the  exception  of 
poor  data  at  IS  ms  delay  (not  plotted),  the  quality  of  the  radar  data 
was  much  improved  over  the  first  day.  This  improvement  was  attributed 
to  the  earlier  starting  time  on  the  second  day,  which  yielded  better 
ionospheric  propagation.  Since  the  radar  data  on  the  second  day  were 
recorded  between  the  hours  of  1620  GMT  and  2003  GMT,  the  wave  directions 
plotted  should  be  applicable  for  mean  winds  blowing  between  the  hours  of 
about  1550  GMT  and  1930  GKT,  respectively. 
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FIGURE  9 ESSA  VIII  CLOUD  PHOTOGRAPH  OVE 
Note  the  cold  front  on  the  left  and  the 


The  magnitude  of  the  wind  was  added  to  the  direction  vectors  in  the 
area  outlined  on  Fig.  9 in  the  vicinity  of  the  cold  front.  The  method 
described  by  Stewart  and  Barnum  [7],  with  the  improvements  illustrated 
by  Maresca  and  Barnum  (8],  was  used.  The  results  are  illustrated  in 
Fig.  10. 


FIGURE  10  WIND  SPEED,  IN  METERS  PER  SECOND,  AND 
DIRECTION  ARROWS  IN  THE  VICINITY  OF 
THE  ATMOSPHERIC  FRONT  ON  26  SEPTEMBER 
1973 


VII  CONCLUSIONS 

The  most  important  contribution  found  in  this  experiment  was  the 
ability  to  detect  the  abrupt  wind  shifts  across  the  cold  front.  Thus, 
when  rapid  spatial  variations  in  wind  fields  occur,  the  high-frequency 
sea  waves  will  also  respond  rapidly  and  this  response  can  be  detected 
by  the  radio  scatter.  High  spatial  resolution  was  thus  essential  for 
the  detection  of  these  abrupt  wind  direction  changes.  In  one  example 
illustrated  herein,  a 16°  shift  occurred  over  only  a half  degree  of  radar 
azimuth,  equivalent  to  about  1_  km  in  cross  range.  Range  resolutions  of 
comparable  size  are  essential  to  the  analysis  of  atmospheric  fronts  and 
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the  gradient  of  wind  stress  across  intense  cyclones.  Resolutions  of 
3 km  are  readily  achievable  with  the  WARF  radar. 

Over  907.  of  high-resolution  sea  backscatter  Doppler  spectra  recorded 
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on  these  two  successive  days  over  an  area  exceeding  3 X 10  (km)  were 
readily  usable  for  calculation  of  the  mean  oceanic  wind  direction.  Due 
to  ionospheric  perturbation  of  the  echoes,  a lesser  number  could  be  used 
for  calculation  of  wind  speed. 

Ionospheric  conditions  varied  both  in  space  and  time.  These  varia- 
tions were  not  abnormal.  Multipath  was  most  severe  at  ranges  exceeding 
3000  km.  This  limitation  was  reduced  by  recording  data  at  the  longer 
ranges  earlier  in  the  day.  It  is  clear  that  the  percentage  of  data 
usable  for  wind  spaed  can  be  increased  by  matching  radar  parameters  to 
ionospheric  skip  parameters  at  each  observation  point.  For  example,  a 
north-south  azimuth  scan  at  a constant  radio  frequency  should  incorporate 
a variable  time-delay  (range)  gate  to  match  north-south  tilts  in  iono- 
spheric electron  density. 

A variety  of  applications  are  evident  for  the  remote  sensing  of  wind 
fields  by  0TH  radar.  The  capability  for  measuring  the  mean  wind  field 
combined  with  a high-resolution  capability  for  detecting  abrupt  shifts 
in  consistent  winds,  should  lead  to  a significant  capability  for  tracking 
and  surveying  hurricanes. 

Future  work  will  be  centered  on  radar  data  analysis  with  good  ground 
truth  under  a variety  of  conditions,  and  on  the  development  of  improved 
processing  and  operating  techniques  in  light  of  ionospheric  variations. 
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Appendix  B 


MEASUREMENT  OF  OCEANIC  WIND  SPEED 
FROM  HF  SEA  SCATTER  BY  SKYWAVE  RADAR 


ABSTRACT 

Remote  measurements  of  the  spatial  mean  ocean  wind  speeds 
were  obtained  using  Doppler  spectra  resolved  to  0.08  Hz  from 
high-resolution  HF  skywave-radar  backscatter  measurements  of  the 
ocean  surface.  A standard  deviation  of  2.4  m/s  resulted  from 
the  correlation  of  observed  winds  over  the  ocean  and  the  broaden- 
ing of  the  Doppler  spectra  in  the  vicinity  of  the  higher  first- 
order  Bragg  line.  This  broadening,  for  Doppler  spectra  unperturbed 
by  the  ionospheric  propagation,  is  proportional  to  the  increase 
in  power  caused  by  higher-order  hydrodynamic  and  electromagnetic 
effects  in  the  vicinity  of  the  Bragg  line  and  inversely  propor- 
tional to  the  square  root  of  the  radio  frequency.  A lower  bound 
on  the  measure  of  wind  speed  was  established  at  5 m/s  by  the 
low-resolution  spectral  processing  and  low  second-order  power. 

An  upper  limit  is  suggested  by  the  steep  slope  in  the  region  of 
the  sea  backscatter  spectrum  outside  the  square  root  of  two 
times  the  first-order  Bragg  line. 


I INTRODUCTION 

High-resolution  HF  OTH-B  radar  measurement  of  ocean  surface  winds 
over  large  areas  of  ocean  are  important  to  the  study  of  large-scale  air/ 
sea  energy  exchanges.  Empirical  estimates  of  the  ocean  surface  wind 
are  an  indirect  measurement  based  on  the  ocean  wave  field  which  acts 
as  a scattering  surface.  Theoretical  models  of  the  hydrodynamic  and 
electromagnetic  contributions  to  the  backscatter  Doppler  spectrum  and 
their  comparison  with  data  obtained  from  HF  surface-wave  radar  have 
been  published  by  Barrick  [1],  and  Johnstone  [2].  The  surface  roughness 
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is  represented  by  an  assumed  form  of  the  ocean  wave  spectrum.  If  the 
assumed  functional  form  of  the  wave  spectrum  results  in  a good  fit 
between  the  model  of  the  Doppler  spectrum  and  the  measured  Doppler 
spectrum,  then  sea  state  is  predictable.  For  fully  developed  condi- 
tions, the  match  between  the  theory  and  the  measured  return  is  good. 
Using  the  second-order  theory,  Johnstone  correlated  wind  speed  with 
the  ratio  of  the  second-order  power  to  the  first-order  power.  Simpler 
determinations  of  wind  speed  from  empirical  correlations  based  on  the 
second-order  theory  have  been  developed  by  Ahearn  et  al.  [3]  and  Stewart 
and  Barnum  [4].  Ahearn  e_t  al..  correlated  the  increase  in  power  measured 
at  the  lowest  level  between  the  Bragg  lines  with  the  wind  speed,  and 
Stewart  and  Barnum  correlated  the  increase  in  the  -10  dB  spectral  width 
of  the  Doppler  spectrum,  defined  as  B,  in  the  vicinity  of  the  higher 
Bragg  line. 

The  objectives  of  this  paper  are  to  report  an  improvement  in  the 
correlation  of  B with  wind  speed,  originally  presented  by  Stewart  and 
Barnum,  by  accounting  for  a difference  in  radar  operating  trequencies , 
and  to  discuss  the  accuracy  of  the  technique. 

II  THEORETICAL  CONSIDERATIONS  AND  DATA  NORMALIZATION 

The  wind  speed  obtained  from  the  spectral  broadening  in  the  vicinity 
of  the  stronger  Bragg  line  is  a spatial  average  of  the  mean  wind  required 
to  generate  the  waves  in  a given  area  of  ocean  at  the  time  of  the  radar 
backscatter  measurement;  it  is  not  a measure  of  the  instantaneous  wind. 

B is  a measure  of  this  spectral  broadening  on  data  processed  with  12.8  s 
of  coherent  integration  with  cosine-squared  weighting  on  the  received 
waveform.  It  was  defined  [4]  as  the  width  of  the  Doppler  spectrum 
10  dB  down  from  the  maximum  Bragg  line  power.  At  this  spectral  resolu- 
tion (0.078  Hz),  the  higher-order  contributions  are  separate  [Fig.  1(a)], 
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SPECTRUM  AMPLITUDE 


DOPPLER  — H» 


FIGURE  1 EXAMPLE  OF  DOPPLER  SPECTRA  COHERENTLY  INTEGRATED  FOR 
(a)  12.8  $ AND  <b)  102.4  s FOR  A GIVEN  RADAR  DWELL 


but  not  uniquely  resolvable  from  the  Bragg-line  power  [Fig.  1(b)].  The 
spectral  broadening  in  this  region  has  been  shown  [1],[2]  to  be  a function 
of  the  higher-order  electromagnetic  and  hydromagnetic  effects  of  scatter- 
ing from  the  ocean  surface.  It  is  probable  that  the  increase  in  second- 
order  power  in  the  Doppler  spectrum  is  a result  of  the  low-frequency 
portion  of  the  wave  spectrum.  This  is  suggested  by  a comparison  of  the 
second-order  power  obtained  from  the  form  of  the  ocean  wave  spectrum 
suggested  by  Phillips  [5],  which  has  a vertical  low-frequency  cutoff 
that  eliminates  the  energy  contribution  at  the  lower  frequencies,  with 
the  form  of  the  ocean  wave  spectrum  suggested  by  Pierson  and  Moscowitz  [6], 
which  includes  the  energy  at  these  lower  frequencies. 

The  resonant  Bragg  scatter  amplitude  from  particular  ocean  waves 

(Doppler  = Af^,  where  n = 1,  2,  3,  4 ...)  saturates  at  particular 

wave  heights  and  no  further  increase  in  power  is  experienced.  For 

example,  sustained  winds  greater  than  12  m/s  will  saturate  the  scatter 

for  n = 1,  2,  and  3 above  a frequency  of  10  mHz.  However,  second-order 

contributions  in  the  region  between  zero  Doppler  and  the  Bragg  line 

3/4 

(Af^)  and  between  the  Bragg  line  and  the  2 Af^  increase  with  increas- 
ine  surface  roughness  and  are  particularly  sensitive  to  increasing  wind 
speeds  [2],  The  theory  shows  that  the  increase  in  power  between  Af^  and 
Af^  is  a function  of  wind  direction  and  wind  speed.  The  greatest 
increase  in  second-order  power  is  found  close  to  the  Bragg  line  as  the 
wind  speed  increases.  For  low  wind  speeds,  less  than  about  5 m/s, 
second-order  contributions  near  the  Bragg  line  are  negligible.  For 
high  wind  speeds,  the  second-order  power  contributions  continue  to 
increase  in  the  region  between  the  Bragg  line  and  Jl  Af^*  The  second- 
order  theory  shows  that  contributions  between  the  positive  and  negative 
first-order  Bragg  lines  are  more  a function  of  wind  direction  and  do 
not  increase  significantly  with  wind  speed  at  a given  direction.  The 
scatter  at  Dopplers  higher  than  the  second-order  contribution  at 
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2 Af  is  not  sensitive  to  an  increase  in  surface  roughness  [2]  and  is 
therefore  not  a good  measure  of  sea  state  and  therefore  wind  speed. 

The  Doppler  frequencies  of  the  Bragg  line  and  surrounding  second- 

order  singularities  are  inversely  proportional  to  the  square  root  of  the 

radio  frequency.  Thus,  the  Doppler  width  of  second-order  scatter  between 
3/4 

Af^  and  2 Af  ^ will  exhibit  a functional  dependence  that  is  also  pro- 
portional to  the  inverse  of  the  square  root  of  frequency.  Also,  as 
evident  in  Barrick  e£  al . [7],  the  magnitude  of  the  second-order  scatter 
relative  to  the  first-order  Bragg  line  will  increase  with  increase  in 
the  radio  frequency  at  a constant  wind  speed.  The  relative  position  of 
the  second-order  peak  also  changes;  thus,  it  is  not  yet  clear  how  the 
overall  frequency  dependence  is  affected. 

The  Bragg  line  can  be  normalized  to  unity  by  dividing  through  the 
Doppler  spectrum  by 

JT* 

V o rtc 

where  f is  the  radio  frequency,  g is  the  acceleration  due  to  gravity, 

and  c is  the  speed  of  light.  Since  g/rc  is  a constant  and  since  the 

mean  frequency  of  the  data  is  approximately  20  MHz,  the  data  were 

1/2 

normalized  by  (20/f  ) , where  f is  the  radio  frequency  in  MHz. 

o o 

III  DATA  ANALYSIS  AND  RESULTS 

The  HF  OTH-B  measurements  were  taken  as  part  of  both  the  NORPAX 
experiment  [4]  and  other  experiments  using  the  Wide  Aperture  Research 
Facility  described  by  Barnum  [8]  and  Barnum  e£  al.  [9],  The  Doppler 
spectrum  used  in  the  B measurement  was  obtained  by  incoherently  averaging 
40  Doppler  spectra  over  a scattering  patch  of  15  km  by  60  lan  for  12.8  s 
of  coherent  integration.  Five  adjacent  areas  of  ocean  (Fig.  2)  separated 
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by  one-quarter  degree  were  simultaneously  processed  in  real  time.  Only 
spectrally  pure  ionospheric  propagation  data  were  analyzed  and  only  the 
minimum  B out  of  the  five  spectra  was  selected  to  use  in  the  correlation. 

The  observed  wind  speeds  were  obtained  from  Research  Vessel  Flip, 

U.S.  Navy  ships,  and  ships  of  opportunity;  therefore,  the  accuracy  of 
the  values  of  wind  speed  observed  in  situ  will  differ.  The  accuracy 
may  range  from  0.5  m/s  to  over  2.0  m/s,  and  depends  on  the  accurate 
determination  of  the  ship's  speed,  the  precision  and  calibration  of  the 
anemometer,  and  the  length  of  the  wind  record  used  to  obtain  an  average 
value . 

The  data  set  is  presented  in  Table  1 and  includes  a subset  of  the 
data  in  [4],  To  eliminate  bias  in  the  correlation,  only  one  B at  each 
geographical  location  was  selected  for  an  observed  wind  on  any  given  day, 
and  no  more  than  three  points  for  any  given  storm. 

A least-squares  fit  was  made  to  the  scatter  diagram  of  the  unnormalized 
B and  observed  wind  speed  for  all  radar  frequencies.  The  resulting  standard 
deviation  was  3.5  m/s,  similar  to  [4].  As  noted  above,  much  of  the  standard 
deviation  can  be  accounted  for  by  the  differences  in  radar  frequency 
associated  with  each  R value. 

1/2 

The  data  were  normalized  by  (20/f  ) , and  are  plotted  in  Fig.  3. 

o 

The  least-squares  fit  to  these  points  is 

U = 33B ' - 2 (1) 

where 

B'  = (20/f  )1/2  B 
o 

and  U is  the  wind  speed  in  m/s.  The  correlation  coefficient  is  0.8, 
and  the  standard  deviation  is  2.4  m/s.  The  previous  influence  of  the 
higher  wind  speeds  on  the  regression  fit  is  not  significant.  The 
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Table  1 


RAW  DATA  USED  IN  CORRELATION  OF  B AND  WIND  SPEED 
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FIGURE  3 CORRELATION  OF  WIND  SPEED  WITH  SPECTRUM  BROADENING 


accuracy  of  the  estimate  of  wind  speed  using  Eq.  (1)  is  on  the  same  order 
of  magnitude  as  the  ship  of  opportunity  reports  presently  used  to  obtain 
the  wind  field  over  a large  area  of  ocean. 


IV  DISCUSSION  OF  ACCURACY 

The  accuracy  of  this  correlation  is  dependent  on  the  natural 
variability  of  the  instantaneous  wind  field,  the  accuracy  of  the 
winds  measured  in  situ,  the  fetch  and  duration  of  sustained  winds, 
and  various  forms  of  spectral  broadening  in  the  vicinity  of  the 
first-order  Bragg  line  that  may  not  be  related  to  the  local  winds 
at  the  scattering  area. 
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A.  Natural  Wind  Variability 

The  standard  deviation  of  the  actual  instantaneous  wind  speed  will 
vary  from  region  to  region,  but  is  usually  assumed  greater  in  the  tem- 
perate zones  than  in  the  tropical  zones  where  the  trude  winds  are  approxi- 
mately constant  for  long  time  periods.  If  the  wind  is  not  uniform  over 
the  scattering  area,  the  observed  wind,  a point-source  measurement, 
will  not  be  representative  of  the  spatially  averaged  wind  speed  obtained 
by  radar  measurements.  In  addition,  the  instrumentation  used  to  measure 
the  winds  in  situ  will  possess  its  own  source  of  error.  Measurements 
from  a moving  platform  will  be  less  accurate  than  from  a stationary  plat- 
form. Errors  may  be  as  great  as  2 m/s.  The  standard  deviation  of  the 
observed  wind  field  determined  from  measurement  error  and  the  natural 
wind  field  variability  define  the  minimum  acceptable  accuracy  of  our 
correlation. 

B.  Fetch  and  Duration 

The  inference  of  the  mean  wind  field  by  ocean  wave  scatter  using  B 
is  dependent  on  fetch  and  duration  [4]  as  well  as  wind  speed.  Since  it 
is  assumed  that  the  mean  wind  field  to  generate  the  ocean  surface  is  the 
mean  wind  averaged  over  that  same  time  period,  only  a fully  developed 
sea  will  give  a unique  result.  The  data  are  insufficient  to  determine 
if  a 20  m/s  wind  was  blowing  for  1 hour  or  if  a 8 m/s  wind  was  blowing 
for  17  hours,  over  a constant  fetch  at  the  time  of  the  backscatter 
measurement.  It  appears  that  the  measure  should  be  strictly  valid  only 
for  a fully  developed  sea.  Since  wave  growth  is  not  considered  in  the 
correlation,  B should  continue  to  increase  until  the  sea  is  fully  de- 
veloped or  at  its  maximum  energy  state  for  the  wind  system.  If  too 
few  values  of  B are  obtained  in  a given  area,  one  cannot  ensure  that  B 
has  reached  its  maximum  value.  This  should  be  responsible  for  some  of 


46 


the  scatter  In  Fig.  3.  We  believe  that  several  values  of  B should  be 
measured  over  time  on  the  order  of  hours,  and  that  the  maximum  value 
of  B,  coincident  with  clean  propagation,  should  be  selected.  It  could 
then  be  assumed  that  the  measured  B has  approached  the  B of  a fully 
developed  sea. 

C.  Unwanted  Spectral  Broadening 

Remaining  potential  errors  result  from  unwanted  spectral  broadening 
in  the  vicinity  of  the  higher  first-order  Bragg  line.  These  include 
ionospheric  multipath  and  fine-scale  motion,  the  possible  influence  of 
swell  and  diverging  (converging)  currents,  and  the  degree  of  statistical 
fluctuation  of  the  spectral  components  in  the  sea  scatter. 

All  OTH-B  data  are  subject  to  ionospheric  perturbation,  and  only 
backscatter  data  free  of  the  effects  of  multipath  and  fine-scale 
ionospheric  motion  in  the  region  of  the  B measurement  can  be  meaningfully 
analyzed.  Five  common  ionospheric  perturbations  to  the  Doppler  spectrum 
are  shown  in  Fig.  4:  (a)  insufficient  signal  strength;  (b)  multipath 

outside  the  B region;  (c)  multipath  inside  the  B region;  (d)  multipath 
in  the  B region  not  easily  distinguishable;  and  (e)  fine-scale  motion 
in  the  B region  (not  shown).  If  the  multipath  does  not  interfere  with 
the  measurement  of  B (Fig.  3(b)],  then  the  measurement  can  still  be 
achieved.  Only  the  return  where  obvious  broadening  has  occurred  but  no 
multipath  is  visible  presents  difficulty  (Fig.  3(d)],  The  five-azimuth 
display  aids  in  assessment.  In  this  case,  the  Doppler  spectra  on  adjacent 
azimuths  have  obvious  multipath  [Fig.  3(d)]  and  it  must  be  concluded  that 
the  center  spectrum  is  contaminated  and  unsuitable  for  analysis.  This 
problem  is  minimal  in  the  correlation  leading  to  (1),  since  only  clean 
spectra  were  used.  To  minimize  ionospheric  broadening  shown  in  Fig.  4(e), 
only  the  minimum  value  of  B was  selected  for  analysis  from  many  Doppler 
spectra  taken  in  a region  at  each  time  period. 
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FIGURE  4 EXAMPLES  OF  DOPPLER  SPECTRA,  (a)  Insufficient  signal  strength,  (b)  Multipath 
(two  paths)  that  does  not  affect  the  measurement  of  B.  (c)  Multipath  (two  paths) 
that  affects  the  measurement  of  B (Hz).  With  further  processing,  the  effects  of  the 
unwanted  path  might  be  subtracted  to  obtain  a clean  Doppler  spectra,  (d)  Multipath 
in  the  adjacent  Doppler  spectra,  simultaneously  recorded  on  adjacent  azimuths, 
indicating  that  the  broadening  in  the  center  Doppler  spectrum  is  due  to  ionospheric 
motion. 
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Spectral  broadeining  of  the  Bragg  line  might  also  occur  from  swell 
[4],  and  from  diverging  or  converging  surface  currents  within  the 
scattering  patch.  Coherent  scatter  from  swell  [2]  appears  at  Dopplers 
well  removed  from  the  Bragg  line  at  10  MHz  (and  higher),  and  therefore 
such  scatter  would  not  broaden  the  spectra.  Modulations  of  the  first- 
order  scatter  by  large  swell  was  not  considered  a source  of  error  in 

[4],  The  narrow  beam  (1/3°)  of  the  WARF  radar  minimizes  the  possibility 

of  diverging  or  converging  currents  occurring  in  the  scattering  area; 
however,  such  currents  could  produce  serious  misinterpretation  of  the 
data  for  larger  antenna  beamwidths. 

To  eliminate  statistical  variation  in  the  spectr  due  to  random  sea 
scatter  amplitudes,  the  experimental  data  were  analyzed  to  determine  the 
optimal  tradeoff  among  spatial  averaging,  and  coherent  and  incoherent 
spectrum  integration.  The  ratio  of  the  approach  and  recede  first-order 
Bragg  lines  defined  as  X [9]  was  used  as  the  parameter  of  the  test. 

Using  a data  record  of  up  to  256  s,  the  analysis  proceeded  for  102,  51, 

26,  and  13  s of  coherent  integration,  and  1,  2,  4,  8,  12,  16,  and  20 
incoherent  averages  for  various  coherent  integration  times.  It  was 
found  that  longer  coherent  integration  time  Increased  spectral  resolu- 
tion but  did  not  reduce  statistical  scatter.  Incoherent  averaging  of 
the  spectra  reduced  the  scatter  (as  would  be  expected  for  this  assumed 
Gaussian  process).  Based  on  the  results  of  this  analysis  (Fig.  5),  it 
was  concluded  that  four  to  eight  incoherent  averages  are  required  per 
range  line  regardless  of  the  coherent  integration  time.  The  incoherent 
spatial  average  of  the  backscatter  measurements  must  be  increased  for 
decreased  incoherent  averaging. 

It  is  concluded  that  the  scatter  in  the  data  reported  here  is  due 
primarily  to  the  natural  variability  of  the  mean  wind  field,  and  the 
error  in  the  in  situ  measurement  of  the  observed  winds  (particularly 
for  the  higher  wind  values). 
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NUMBER  OF  INCOHERENT  AVERAGES  — map* 


FIGURE  5 THE  EFFECT  OF  INCOHERENT  AVERAGING  ON  OBTAINING  A MEAN 
VALUE  FOR  X.  Doppler  spectra  coherently  integrated  for  12.8,  25.6, 
51.2,  and  102.4  s were  averaged  to  obtain  a mean  value  of  X.  The  mean 
value  of  X M dB  is  approached  for  4 to  8 incoherent  averages  (maps). 


V RANGE  OF  B 

A minimum  value  of  B is  determined  by  the  processing.  According  to 
our  present  automatic  measure  of  B with  a 12.8-s  integration,  its  minimum 
value  for  no  second-order  contribution  ranges  from  0.197  Hz  to  0.225  Hz, 
depending  on  the  location  of  the  first-order  Bragg  line  relative  to  the 
nearest  Doppler  cell  generated  by  the  FFT.  These  minimum  values  of  B 
above  were  calculated  for  the  Bragg  line  centered  on  a Doppler  cell  and 
located  midway  between  Doppler  cells,  respectively.  A B of  0.2  Hz  cor- 
responds to  a wind  speed  of  5 m/s  at  20  MHz  (Fig.  3).  Choosing  the  minimum 
value  of  B from  a number  of  relatively  closely  spaced  samples  following 
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unperturbed  propagation  tends  to  result  in  selection  of  Doppler  spectra 
in  which  the  Bragg  line  is  centered  or  is  close  to  centered  on  a Doppler 
cell . 

At  the  other  extreme,  the  theoretical  spectra  calculated  by  Johnstone 
[2]  suggest  an  upper  limit  on  the  wind-speed  calculation  using  B.  As  the 
wind  speed  increases,  the  power  in  the  spectral  peak  found  at  Jl  Af  ^ 
also  increases.  If  this  peak  saturates  below  -10  dB  the  linear  approxi- 
mation on  which  the  correlation  is  based  is  reasonable.  When  this  peak 
saturates  above  -10  dB,  B increases  only  very  slowly  with  increasing  wind, 
due  to  the  steep  slope  of  the  region  outside  this  spectral  component,  and 
the  linear  approximation  to  U(B)  will  become  Invalid.  Although  B will 
continue  to  increase  with  increasing  wind  speed,  this  Increase  will  not 
be  detectable  by  this  method. 


VI  CONCLUSION 

The  spectral  broadening  in  the  vicinity  of  the  first-order  Bragg 

line  at  a constant  sea  state  is  a function  of  the  radio  frequency.  One 

factor  of  this  function  is  inversely  proportional  to  the  square  root  of 

-1/2 

frequency.  Normalizing  the  data  by  f resulted  in  a standard  devia- 
tion of  2.4  m/s,  which  is  equal  to  the  ship  of  opportunity  data.  The 
correlation  is  valid  for  wind  speeds  greater  than  5 m/s  at  20  MHz  (Fig.  3); 
an  upper  wind-speed  limit  is  suggested  by  theory,  but  is  unknown.  The 
lower  limit  is  determined  by  the  minimum  spectral  width  possible  in  the 
spectral  processing  and  the  negligible  second-order  contributions  near 
the  first-order  Bragg  line  for  winds  less  than  5 m/s.  The  use  of  B to 
infer  oceanic  winds  is  particularly  suited  to  covering  large  ocean  areas 
in  detail,  since  only  12.8  s of  coherent  integration  are  required  for 
the  measurement.  It  is  especially  noteworthy  that  the  relatively  low 
wind-speed  error  (2.4  m/s)  was  found  using  B,  even  though  winds 
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from  5 to  25  m/s  were  observed,  and  the  fetch  and  duration  of  the  winds 
were  unknown.  Further  work  will  confirm  whether  or  not  this  was  pure 
coincidence. 
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Appendix  C 


MEASUREMENT  OF  AN  OCEAN  WAVE  SPECTRUM 
BY  HF  SKYWAVE  RADAR 


ABSTRACT 

An  ocean  wave  spectrum  for  the  North  Facific  Ocean  (29°N, 
137 °W)  was  inferred  from  sea  backscatter  measurements  taken 
by  HF  skywave  radar  on  12  February  1975.  A trial-and-error 
solution  was  used  to  estimate  this  spectrum  by  comparing  a 
theoretical  second-order  radio-wave  scattering  model  of  the 
Doppler  spectrum,  produced  by  assuming  an  ocean  wave  spectrum, 
to  the  observed  Doppler  spectrum  until  a match  was  found. 

The  comparison  of  the  final  calculated  Doppler  spectrum  and 
the  measured  Doppler  spectrum  showed  fair  agreement.  The 
significant  wave  height  derived  using  this  technique  agreed 
to  within  ten  percent  of  the  in  situ  observations  of  the  wave 
height. 


I INTRODUCTION 

High-frequency  groundwave  radar  measurements  have  demonstrated  the 
feasibility  of  determining  the  fully  developed  ocean  wave  spectrum  from 
sea  backscatter.  The  rangj  of  groundwave  radar  measurements  is  limited 
to  less  than  200  km.  Eventual  application  of  HF  measurements  will  be 
at  skywave  radar  ranges  extending  out  to  3000  km. 

The  second-order  contributions  to  the  Doppler  spectrum  obtained  from 
HF  sea  backscatter  measurements  were  modeled  theoretically  by  Barrick  [1] 
and  Johnstone  [2].  Assuming  three  forms  of  the  fully  developed  ocean 
wave  spectrum  suggested  by  Phillips  [3],  Dierson  and  Moscowitz  [4],  and 
Tyler  et  al . [5],  Johnstone  found  fair  agreement  between  the  theoretical 

Preceding  page  blank 
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Doppler  spectra  and  the  observed  Doppler  spectra  obtained  from  the  ground 
wave  radar  experiments  conducted  at  Wake  Island  by  Tyler  et  al.  [5]. 

The  best  agreement  was  found  for  the  Pierson-Moscowitz  form  of  the  wave 
spectrum. 

It  is  essential  that  an  ocean  wave  spectrum  can  be  determined 
uniquely  from  the  backscatter  data  without  assuming  the  sea  is  fully 
developed.  A trial -and-error  solution  was  employed  to  match  the 
theoretical  Doppler  spectrum,  based  o.i  an  assumed  ocean  wave  spectrum, 
to  the  observed  Doppler  spectrum.  The  purpose  of  this  paper  is  to 
demonstrate  the  potential  of  using  an  HF  skywave  radar  to  remotely 
measure  the  ocean  wave  spectrum  at  distances  up  to  3000  km. 


II  THEORY 

The  theoretical  model  of  second-order  HF  radio-wave  scattering  from 
the  sea  derived  by  Johnstone  [2]  to  determine  the  incremental  surface 
radar  cross  section  per  unit  area  is  used  to  estimate  the  ocean  wave 
spectrum  from  HF  skywave  backscatter  measurements.  The  radar  cross 
section  per  unit  area,  a,  is 


c 


(1) 


where  A is  the  surface  area  of  illumination,  R is  the  distance  of  the 
transmitter  and  receiver  from  the  target,  E is  the  scattered  electric 
field  measured  at  R,  and  E^  is  the  electric  field  incident  at  the  ocean 
patch.  The  effect  of  the  scattering  surface  on  the  received  power  is 
expressed  by  the  radar  cross  section.  The  radar  cross  section  per  unit 
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area  per  unit  frequency  for  bistatic  radar  and  vertical  polarization 
expanded  to  first,  second,  and  higher  order  is 

a (id)  ■ (to ) + (ou)  + higher  order  terms  (2) 

w w w 

where  id  is  the  radian  radio-wave  frequency.  The  derivation  of  the  ex- 
pressions for  Eq.  (2)  is  found  in  Johnstone  [2].  Proportional  expressions 
that  show  the  relationship  of  the  ocean  wave  spectrum  to  the  radar  cross 

section  per  unit  area  per  unit  frequency  are  given,  (tu),  the  first- 

w 

order  term  or  Bragg  line  response,  represents  a resonant  interaction  of 
the  ocean  waves  of  length  L equal  to  one-half  the  radio  waves  of  length 
\ and  is  proportional  to  the  ocean  wave  spectrum,  W(k). 

o™  («)  « W , 0^  6(T]  ± o)B)  (3) 

where  6(T]  ± u>  ) is  the  Dirac  delta  function,  T|  is  the  scatter  radio-wave 
B 

Doppler  shift,  and  id  is  the  ocean-wave  radian  frequency  for  which  first- 
order  Bragg  scattering  occurs.  Grazing  incidence  and  monostatic  geometry 
are  assumed  for  the  radar. 

The  amplitude  of  the  Bragg  lines  in  the  theoretical  Doppler  spectrum 

calculated  from  Eq.  (3)  is  not  comparable  to  the  amplitude  of  the  Bragg 

lines  in  the  observed  Doppler  spectrum  because  the  former  is  power  while 

(2) 

the  latter  is  power  spectral  density,  a (u>)  is  the  second-order  term 

which  results  in  the  Doppler  continuum  and  is  directly  comparable  to  the 

observed  Doppler  spectrum.  Therefore,  only  the  second-order  contributions 

(2) 

to  the  Doppler  spectrum  are  compared,  (id)  is  related  to  the  ocean 

wave  spectrum 

oo 

’ j*  if  + 4r")  M11  ‘ "l  + V “<"1>  “h  ^ <4> 

o 
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where  I"  and  r are  the  electromagnetic  and  hydrodynamic  kernels, 
EM  H 

respectively,  and  k is  the  radian  ocean  wave  number. 


Ill  THEORETICAL  MODEL  OF  HF  SECOND-ORDER  SCATTERING 
FROM  THE  OCEAN  SURFACE 

The  second-order  radio-wave  scattering  model  described  by  Johnstone 
[2]  assumes  a fully  developed  form  of  the  ocean  wave  spectrum.  Since 
the  sea  may  be  growing  or  decaying  at  the  time  of  the  backscatter  measure- 
ments, the  form  of  the  ocean  wave  spectrum  is  unknown.  Practical  applica- 
tion of  this  model  requires  that  for  a given  measured  Doppler  spectrum, 
a unique  wave  spectrum  can  be  calculated.  The  required  mathematical 
solution  for  this  inversion  operation  has  not  been  developed,  however, 
and  a trial-and  error  solution  is  postulated.  First,  a form  of  the 
ocean  wave  spectrum  is  assumed.  Second,  a theoretical  Doppler  spectrum 
is  calculated  based  on  this  ocean  wave  spectrum.  Third,  the  theoretical 
Doppler  spectrum  in  Eq.  (4)  is  compared  to  the  radar-measured  Doppler 
spectrum.  If  the  agreement  is  good,  then  the  iterative  solution  is 
terminated.  If  not,  adjustments  are  made  in  the  assumed  ocean  wave 
spectrum  and  another  iteration  is  performed. 

The  trial -and-error  solution  progressed  in  three  steps.  First, 
the  fully  developed  form  of  the  ocean  wave  spectrum  suggested  by  Pierson 
and  Moscowitz  [4]  was  used: 


W(w) 


0^78 


(5) 


where  u is  the  wind  speed  and  g is  the  acceleration  due  to  gravity.  The 
theoretical  Doppler  spectra  obtained  using  Eq.  (5)  to  describe  the  surface 
toughness  were  compared  to  the  observed  Doppler  spectra.  Second,  if 
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satisfactory  agreement  was  not  found,  then  the  shape  of  the  ocean  wave 
spectrum  was  changed  for  an  assumed  constant  value  of  the  amplitude  and 
frequency  of  the  spectral  peak.  This  was  accomplished  by  normalizing 
Eq.  (5)  to  obtain 


W(uu  ± ou  ) 

o 


A_ 

Nl 


(6) 


where  A,  B,  Nl,  and  N2  are  constants  that  can  be  varied  to  change  the 

shape  of  the  spectrum,  |wp|  is  the  peak  amplitude  of  the  wave  spectrum 

in  Eq.  (5),  |sp|  is  the  assumed  amplitude  of  the  spectral  peak  of  the 

wave  spectrum,  and  uj^  is  a shift  in  the  ocean  wave  frequency  from  the 

peak  frequency  in  Eq.  (5).  A directional  distribution  for  ocean  gravity 
2 

waves  of  cos  (9/2)  was  assumed. 

The  theoretical  Doppler  spectrum  is  compared  to  the  observed  Doppler 

spectrum  by  matching  the  Bragg-line  frequencies,  Af,  and  the  peak  ampli- 

l/2  3 /4 

tudes  of  the  second-order  contribution  at  2 Af^  and  2 Af ^ . This 

is  necessary,  since  the  amplitudes  of  both  the  theoretical  and  observed 
Doppler  spectra  are  in  relative  units.  The  fit  was  judged  visually, 
according  to  the  criteria  described  below. 

The  significant  wave  height,  was  estimated  from  Eq.  (5)  or 

Eq.  (6)  by 

OD 

— 4a  - j W(uu)  du;  (7) 

0 

where  c is  the  standard  deviation  of  the  wave  record  based  on  a Rayleigh 
distribution  function  [6]. 
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SPECTRUM  AMPLITUDE 


IV  SKYWAVE  SEA  BACKSCATTER  MEASUREMENTS 


High-frequency  sea  backscatter  measurements  relatively  free  of 

ionospheric  perturbations  were  taken  by  WARF,  a skywave  radar  located  in 

central  California  on  12  February  1975  as  part  of  the  NORPAX  experiment 

and  other  experiments.  The  data  were  obtained  from  a spatial  average 

2 

of  the  ocean  waves  in  a 900-km  area  located  in  the  North  Pacific  Ocean 
at  29°N,  137  °W.  The  radio-wave  propagation  was  via  a stable  sporadic-E 
layer  over  the  1600-km  distance  from  the  radar  to  the  illuminated  ocean 
patch . 

The  sea  backscatter  data  were  analyzed  by  a range/Doppler  processing 
technique.  Four  individual  range  resolution  cells  were  processed  for 
102.4  s of  integration  to  obtain  a Doppler  spectrum  for  each  radar  dwell. 
Two  such  Doppler  spectra  were  incoherently  averaged  to  obtain  the  reduced 
data  shown  in  Fig.  1.  In  general,  the  skywave  data  are  excellent,  with 
minimal  ionospheric  perturbations. 


FIGURE  1 EXAMPLE  OF  THE  OBSERVED  HF  SEA  BACKSCATTER  MEASUREMENTS 
TAKEN  IN  THE  NORTH  PACIFIC  OCEAN 
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V OCEAN  WAVE  OBSERVATIONS 


Ocean  wave  spectra  were  not  available  for  this  experiment.  Visual 
wave-height  observations  of  2.1  m and  1.5  m were  obtained  from  two  Navy 
ships  operating  in  the  experimental  area.  The  empirical  wave-hindcastlng 
curves  developed  by  Bretschneider  [6]  were  used  to  predict  the  significant 
wave  height.  Two  estimates  of  the  surface  wind  were  used--corrected 
surface  geostrophic  wind,  and  hourly  wind  reports  obtained  from  the  two 
Navy  ships.  The  results  are  shown  in  Table  1.  The  estimates  by  all 
methods  are  within  10  percent. 


Table  1 

OBSERVED  AND  HINDCAST  OCEAN  WAVE  HEIGHT 


Method 

Corrected 

Surface 

Geostrophic 

Wind 

(m/s) 

Mean 

Observed 

Wind 

(m/s) 

Significant 

Wave 

Height 

(m) 

Observed 

Wave 

Height 

(m) 

9.0 

10.3 

1.8 

1.9 

Visual 

Observation 

1.8 

VI  RESULTS 

Comparisons  of  the  second-order  theoretical  contributions  to  the 
Doppler  spectra  and  observed  Doppler  spectra  were  made  using  the  assumed 
ocean  wave  spectra  given  in  Table  2 and  shown  in  Fig.  2.  These  comparisons 
are  shown  in  Fig.  3.  The  best  fit  between  the  theoretical  and  experimental 
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Doppler  spectra  occurred  for  the  ocean  wave  spectrum  shown  in  Fig.  3(g) 
for  waves  approaching  the  radar  at  160°  from  the  radar  bearing  for 


S(uu  ± w ) 
o 


0.78 


U) 


0.74 

4 

(JU 


(7) 


where  id  = 0.000.  The  visual  fit  was  based  on  the  match  among  the 
o 

theoretical  and  observed  cutoff  values,  slope  of  the  Doppler  continuum, 
3/4 

and  power  from  2 Af^  to  Af  ^ • Since  the  comparison  between  the  theo- 
retical and  observed  Doppler  spectra  is  made  in  dB,  the  maxima  of  the 
spectral  fluctuations  should  be  weighted  more  than  the  minima. 


The  comparison  of  the  theoretical  and  experimental  Doppler  spectra 

3/4 

is  relative  to  the  spectral  peak  at  2 Af.  The  amplitude  of  the 

1/2  3 //) 

observed  peaks  at  2 Af^  and  2 Af^  may  range  ±2  dB,  and  the  ampli- 
tude of  the  theoretical  peaks  are  subject  to  calculation  errors  since 
they  occur  as  discontinuities  in  the  mathematical  solution.  Therefore, 
the  fit  of  the  entire  second-order  Doppler  continuum  was  weighted  more 
heavily  in  determining  the  best  fit.  The  theoretical  Doppler  continuum 
outside  the  lower  Bragg  lines  is  slightly  lower  than  predicted.  A 
significant  wave  height  of  1.9  m was  calculated  for  the  ocean  wave 
spectrum  given  by  Eq.  (7),  and  is  within  10  percent  of  the  observed 
and  hindcast  wave  heights. 


The  uniqueness  of  Eq.  (7)  is  unknown.  It  can  be  argued  from  the 
trial -and-error  calculation  that  other  ocean  wave  spectra  approximately 
equal  in  total  energy  can  be  derived  using  this  technique.  We  attribute 
the  close  fit  of  several  of  the  theoretical  Doppler  spectra  to  the 
observed  Doppler  spectra  to  the  spectral  fluctuations  in  the  observed 
da  ta . 
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For  example,  Che  fit  between  Che  theoretical  Doppler  spectrum 
derived  from  the  ocean  wave  spectrum  labeled  Fig.  2(a)  in  Table  2 and 
the  observed  Doppler  spectrum  is  also  good.  We  did  not  choose  this 
fit  since  the  second-order  contributions  near  the  higher  Bragg  line 
appeared  too  high.  On  the  other  hand  we  can  easily  eliminate  the  fit 
for  ocean  wave  spectra  of  more  energy  than  Fig.  2(a)  [Fig.  3(f)J  and 
less  energy  than  Fig.  2(a)  [Fig.  3(d)J. 


The  theoretical  Doppler  spectrum  is  particularly  sensitive  to  the 
wave  energy  at  each  frequency  as  well  as  the  total  energy.  The  first 
four  theoretical  Doppler  spectra  shown  in  Figs.  3(a)  through  3(d)  were 
generated  by  an  assumed  form  of  the  ocean  wave  spectra  suggested  by 
Pierson  and  Moscowitz  [4]  for  mean  wind  speeds  of  7.8  m/s,  8.7  m/s, 

9.8  m/s,  and  11.8  m/s.  The  theoretical  Doppler  spectra  were  influenced 
by  the  increased  wave  energy  at  the  lower  frequencies.  The  next  two 
theoretical  Doppler  spectra  shown  in  Figs.  3(e)  and  (f)  were  generated 
by  a narrow  peaked  ocean  wave  spectrum  and  a broad  ocean  wave  spectrum, 
respectively,  at  the  same  peak  amplitude  and  frequency  equal  to  the 
Pierson  and  Moscowitz  form  of  the  ocean  wave  spectrum  for  an  11.8-m/s 
mean  wind.  The  differences  in  the  resulting  Doppler  spectra  are  due 
to  the  difference  in  total  wave  energy  at  the  various  wave  frequencies. 
It  is  interesting  to  consider  identlal  ocean  wave  spectra  at  three 
different  peak  wave  frequencies  [Figs.  2(g),  (h),  and  (i)].  An  ocean 
wave  spectrum  shifted  to  lower  frequencies  results  in  a lower  cutoff 
value  and  decreased  energy  in  the  Doppler  continuum  including  the  peaks 

An  ocean  wave  spectrum  shifted  to  the  higher 


1/2  r . 3/4  . 

at  2 Af ^ and  2 Af 


frequencies  results  in  a higher  cutoff  value  and  increased  energy  in 

X / 2 3 /4 

the  Doppler  continuum  including  the  peaks  at  2 Af^  and  s Af^» 

The  slope  and  energy  content  of  the  theoretical  Doppler  spectra  change 
with  contributions  to  the  ocean  wave  spectrum  at  different  frequencies. 
For  a given  shape  of  the  ocean  wave  spectrum  in  which  the  total  wave 
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I 

energy  may  change,  the  theoretical  Doppler  spectrum  increases  in  power 
but  does  not  change  slope. 

VII  CONCLUSIONS 

Long  dwells  necessary  to  obtain  sufficient  resolution  in  the  Doppler 
spectrum  can  be  obtained  by  skywave  measurements.  We  ha/e  demonstrated 
a trial-and-error  solution  method  to  estimate  the  ocean  vave  spectrum. 

A comparison  between  the  significant  wave-height  calculated  from  the 
ocean  wave  spectrum  and  the  observed  and  hindcast  significant  wave 
heights  showed  agreement  to  within  ten  percent.  Future  work  at  deter- 
mining the  sensitivity  of  the  measurement,  the  effects  of  swell  generated 
from  several  distant  storms,  and  the  effects  of  the  directional  properties 
of  ocean  waves  is  warranted. 
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